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Abstract 

We present an analysis of the Belle measured moments of the lepton energy and hadronic mass 
spectra m. B ^ Xciv decays and the photon energy spectrum m. B ^ Xg^ decays using theoretical 
expressions derived in the IS and kinetic schemes. The magnitude of the Cabibbo-Kobayashi- 
Maskawa matrix element Fcfe, the 6-quark mass and other non-perturbative parameters are ex- 
tracted. In the IS scheme analysis we find \Vcb\ = (41.49 ± 0.52(fit) ± 0.20(rB)) x 10"^ and 
ml^ = (4.729 ±0.048) GeV. In the kinetic scheme, we obtain \Vcb\ = (41.93 ± 0.65(fit) ± 0.07(as) ± 
0.63(th)) X 10-=^ and mf"" = (4.564 ± 0.076(fit) ± 0.003(as)) GeV. 
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I. INTRODUCTION 



The most precise determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ment iV^cbl [1] are obtained using combined fits to inclusive B decay distributions [2, 3]. These 
analyses are based on calculations of the semileptonic decay rate in the frameworks of the Op- 
erator Product Expansion (OPE) [4] and the Heavy Quark Effective Theory (HQET) [2, 5] 
which predict this quantity in terms of \ Vcb\, the 6-quark mass m^, and a few non-perturbative 
matrix elements that enter at the order 

Several studies have shown that the spectator model decay rate, in which bound state 
effects are neglected, is the leading term in a well-defined expansion controlled by the param- 
eter Aqcd/''716 [5-8]. Non-perturbative corrections to this leading approximation arise only 
to order I /ml. The key issue in this approach is the abihty to separate non-perturbative 
corrections (expressed as a series in powers of l/m^), and perturbative corrections (expressed 
in powers of a^). There are various different methods to handle the energy scale /i used to 
separate long-distance from short-distance physics. 

The coefficients of the l/m^ power terms are expectation values of operators that include 
non-perturbative physics. In this framework, non-perturbative corrections are parameterised 
by quark masses and matrix elements of higher dimensional operators which are presently 
poorly known. The experimental accuracy already achieved, and that expected from larger 
data sets recorded by the -B-factories, make the ensuing theory uncertainty a major limit- 
ing factor. The extraction of the non-perturbative parameters describing the heavy quark 
masses, kinetic energy of the b quark and the corrections directly from the data has 

therefore become a key issue. 

The non-calculable, non-perturbative quantities are parametrised in terms of expecta- 
tion values of hadronic matrix elements, which can be related to the shape (characterised 
by moments) of inclusive decay spectra [2, 8, 9]. High precision comparison of theory 
and experiment requires a precise determination of the heavy quark masses, as well as the 
nonperturbative matrix elements that enter the expansion. These are Ai,2 or //^^g which 
parameterise the nonperturbative corrections to inclusive observables at 0(AQ(-,j-)/m^). At 
0(AQQj-)/m^), more matrix elements occur, denoted by pi,2 and ti_4 or pd,ls- 

In this paper we make use of the Heavy Quark Expansions (HQEs) that express the 
semileptonic decay width Fs.i., moments of the lepton energy and hadronic mass spectra in 
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B — > Xciv decays and the photon energy spectrum in S — > Xs^ decays in terms of the 
running kinetic quark masses m^'" and m^™ as well as the IS 6-quark mass m^^. Further 
details of these two schemes are discussed in later sections of the paper. These schemes 
should ultimately yield consistent results for \Vcb\- The precision of the 6-quark mass is also 
important for the determination of \Vub\- the least well understood element in the CKM 
matrix, and a limiting factor in the uncertainty on the unitarity triangle. 

II. EXPERIMENTAL INPUT 

Belle has measured the partial branching fractions B{B — > Xciiy)E(>E^i^ and the first, 
second, third and fourth moments of the truncated electron energy spectrum in B ^ Xceu, 
{E,)e,>e^,^, {{Ee-{Ee)f)E,>E^,^, {{E,-{E,)f)E,>E^,^ and {{E,-{E,) Y) e,>e^,^. for nine dif- 
ferent electron energy thresholds {E^iry = 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8 and 2.0 GeV) [10]. 

We have measured the first, second central and second non-central moments of the hadron 
invariant mass squared (M|-) spectrum in S ^ Xclv, {Mx)e^>e^i^, {{^x - {^x)Y)Ee>E^i^ 
and {Mj^)E(>E^i^ for seven different lepton energy thresholds (-Emm = 0.7, 0.9, 1.1, 1.3, 1.5, 
1.7 and 1.9 GeV) [11]. 

For B we have measured the first and second moments of the truncated photon 

energy spectrum , {E^) E^>E^i^ and {{E^ — {E^)y) E^>E^in- These measurements are available 
for six minimum photon energies (-Emin — 1-8, 1.9, 2.0, 2.1, 2.2 and 2.3 GeV) [12]. 

Hence, there are a total of 71 Belle measurements of inclusive spectra available for use 
in the fits [13]. The measurements used in the IS and kinetic scheme fit analyses are shown 
in Table I. We have excluded measurements that do not have corresponding theoretical 
predictions. Measurements with higher cutoff energies (i.e. electron energy and hadron mass 
moments with -Emin > 1-5 GeV and photon energy moments with E'min > 2 GeV) are not 
used to determine the HQE parameters, as theoretical predictions are not considered reliable 
in this region. Finally, we have also excluded points where correlations with neighbouring 
points are too high. 

The value of 1 is dependent on the B meson lifetimes. The measured semileptonic 
ratios can be written as Bs.\. = TcfrFs.i. in terms of an effective lifetime, Tcs = f+T+ + /oTo- 
Using the most recent world average values for the lifetimes and the 6-hadron fractions, we 
obtain t^s — (1.585 ± 0.006) ps [14]. For simplicity we refer to this quantity simply as tb- 
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IS scheme 


kinetic scheme 


Lepton moments {E^)E^i^ 


n = Emin = 0.6, 1.0, 1.4 
n = 1 Emin = 0.6, 0.8, 1.0, 1.2, 1.4 
n = 2 E^i^ = 0.6, 1.0, 1.4 
n = 3 ^^min = 0.8, 1.2 


n = -Bmin = 0.4, 0.8 
n = 1 Emin = 0.4, 0.8, 1.0, 1.2 1.4 
n = 2 E^^ = 0.4, 0.8, 1.0, 1.2 1.4 
n = 3 Eram = 0.4, 0.8, 1.0, 1.2 1.4 


Hadron moments {M]p) e^-,,^ 


n = 1 £;min = 0.7, 1.1, 1.3, 1.5 
n = 2 Eynin = 0.7, 0.9, 1.3 


n = 1 E^i^ = 0.7, 0.9, 1.1, 1.3 
n = 2 E^min = 0.7, 0.9, 1.1, 1.3 


Photon moments e^^,^ 


n = 1 ^^min = 1.8, 2.0 
n = 2 = 1.8, 2.0 


n = 1 -Bmin = 1-8, 1.9, 2.0 
n = 2 E^i^ = 1.8, 1.9, 2.0 



TABLE I: Experimental input used in the IS and kinetic scheme analyses. The values of -Bmin are 
given in GeV. The IS (kinetic) scheme analysis uses a total of 24 (31) measurements. 

III. IS SCHEME ANALYSIS 
A. Theoretical input 

The inclusive spectral moments of -B — > Xciiy decays have been derived in the IS scheme 
up to 0{l/ml) [2]. The theoretical expressions for the truncated moments have the following 
form (where {X)E^i^ represents any of the experimental observables) : 

{X)e^,^ = + X(2)A + X(3)A2 + X^^^A^ + X(5)Ai + X(«)AAi + X^'^X2 + X^^^AX^ + X^'^p, 
+X(io)p2 + ^("Vi + + X(i3)r3 + X(i^)r4 + X^^'^e + X(i6)^2^^ + ^(17)^^ _ 

(1) 

The coefficients X^''\ determined by theory, are functions of E^am- The non-perturbative 

corrections are parametrized by A (0(mf,)), Ai and A2 (0(l/m^)), and ri, T2, T3, T4, pi and 
P2 {0{l/ml)). 

Predictions for the partial branching fractions are obtained using the following expression, 
B(B XMe.... = {XUeJ^^^S^^^ , (2) 

where {X)B,E^i^ is an expression of the form of Eq. 1, m is the IS reference mass, m — 
mr(is)/2, Vqed = 1-007, and G%my {192Tr^) = 5.4 x 10"^^ 
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In this analysis, we determine a total of seven parameters: \Vcb\, A, Ai, ti, T2, T3 and pi. 
One of the higher order parameters, T4 is set to zero , and from available constraints, e.g. 
B* — B mass splitting, the remaining parameters in Eq. 1 are set to: A2 — 0.1227 — 0.0145Ai 
and P2 = 0.1361 + r2, following advice from Ref. [2]. The parameter A is the difference 
between the 6-quark mass and the reference value about which it is expanded, i.e., A = 
"^T(i5)/2 — ml^ . We will present our results in terms of ml^ in place of A. 



B. The function 

The fit takes into account both experimental and theoretical uncertainties. Following the 
approach in Ref. [2] , an element of the combined experimental and theoretical error matrix 
is given by 

(^ij = cri(^jCij, (3) 

where i and j denote the observables and Cij is the experimental correlation matrix element. 
The total error on the observable i is defined as 

ai = ^J{(7T^Y + {Afr^rri^^Y + (5^/2)2 for the nth hadron moment , (4) 
(Ti = ^((jfP)2 + {Afn{mB/2YY + {Bi/2f for the nth lepton moment , (5) 
Gi = + {Afn{mB/2YY + (5^/2)2 for the nth photon moment , (6) 

and fo — fi — 1, f2 — 1/4: and /a = l/(6-\/3). Here, af^"^ are the experimental errors, Bi — 
X^^^^ are the coefficients of the last computed terms in the perturbation series (providing the 
error on the uncalculated higher order perturbative terms) . The dimensionless parameter A 
contains various theoretical errors (uncalculated power corrections, uncalculated effects of 
order (Q;s/47r)Aqcj-,/m^, and effects not included in the OPE, i.e., duality violation), and is 
multiphed by dimensionful quantities. We fix ^4 = 0.001 as in Ref. [2]. For B Xg^, the 
accessible phase space is limited, and the theoretical extraction of is affected by shape 
function effects. So, A is multiplied by the ratio of the difference from the end point relative 
to E^i^ = 1.8 GeV. 

As the fit does not provide strong constraints on the C(l/mf) parameters, it is necessary 
to provide constraints to ensure their convergence to sensible values. We achieve this by 
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introducing extra terms in the function of the fit, 

{0 \0\<ml, 

where {m^,M^) are both quantities of order Aqcd, and (O) represents any 0{l/ml) param- 
eters. In the fit, we take M^. — — 500 MeV after Ref. [2]. The parameter may have 
a value anywhere between 500 MeV and 1 GeV. 

The overall form of the function used in the IS fit is 

x' = E((^)r"- - • - {x)f) +i2xl+i2 4 , (8) 

i,j i=l i=l 

where {X)f^^^- are the measured moments and (X)]^ are the corresponding IS scheme 
predictions. 

C. Fit results and discussion 

Minimizing the function in Eq. 8 using MINUIT [15], we find the following results for 
the fit parameters, 

\Veb\ = (41.49 ± 0.52fit ± 0.20^^) X 10-^ , 
= (4.729 ± 0.048) GeV , and 
Ai = (-0.30 ± 0.04) GeV^ . 

The first error is the uncertainty from the fit including experimental and theory errors, and 
the second error (on \Vcb\ only) is due to the uncertainty on the average B lifetime. The 
correlations between these fit parameters are provided in Table II. Using the measurement of 
the partial branching fraction at £'min = 0.6 GeV, we obtain for the semileptonic branching 
ratio (over the full lepton energy range), 

B(B XJtv) = (10.62 ± 0.25)% . 

The measured moments compared to the IS scheme predictions are shown in Figs. 1 and 2. 

We assess the stabihty of the fit in two ways (Table III): by repeating the fit only to 
B — > Xc^v data (21 measurements), (a) to (c); and by releasing the constraint on 
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\Vcb\ mf A 



\Vcb\ 1.000 -0.539 -0.330 



mf 1.000 0.871 



1.000 



TABLE II: Correlation coefficients of the parameters in the IS fit. 
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FIG. 1: Fit results to the electron energy spectrum moments. The yellow band represents the 
fit error, and the red band gives the theory and fit errors combined. Filled circles represent data 
points used in the fit, and open circles are points not used in the fit. 
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FIG. 2: Fit results to the hadron invariant mass squared and photon energy spectrum moments. 
The yeUow band represents the fit error, and the red band gives the theory and fit errors combined. 
Filled circles represent data points used in the fit, and open circles are points not used in the fit. 



the higher order parameters, (c) and (f). To study the effect of the estimated theoretical 
uncertainties, we repeat the fit with all theoretical uncertainties set to zero; (b) and (e). All 
studies give consistent results with acceptable values of x^/ndf. Figure 3 shows the contour 
plots for the fits corresponding to Table Ill(a) {B X^^z/ data only) and Ill(d) (full fit). 
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Data used 






2 

'^theory 


xVndf. 


\Vch\ X 103 


mf (GeV) 


Ai (GeV2) 


(a) -£^e(0,l,2,3), 


^1(1,2) 




0.5 


yes 


4.9/13 


41.52 ±0.77 


4.723 ±0.103 


-0.307 ±0.083 


(b) -E^e(0,l,2,3)> 


^1(1,2) 




0.5 


no 


9.4/13 


41.64 ±0.58 


4.655 ± 0.075 


-0.348 ±0.052 


(c) ^e(0,l,2,3)) 


^1(1,2) 




0.8 


yes 


1.9/13 


42.82 ± 1.09 


4.588 ±0.118 


-0.505 ±0.271 


(d) -£^e(0,l,2,3)) 


^7(1,2)' 


^1(1,2) 


0.5 


yes 


5.7/17 


41.49 ±0.52 


4.729 ± 0.048 


-0.302 ± 0.043 


(e) ^e(0,l,2,3)> 


-£'7(1,2)) 


^1(1,2) 


0.5 


no 


10.8/17 41.42 ±0.45 4.695 ± 0.046 


-0.321 ±0.035 


(f) -£^6(0,1, 2,3): 


-^7(1,2), ^1(1,2) 


0.8 


yes 


3.9/17 


42.19 ±0.73 


4.709 ± 0.066 


-0.489 ± 0.087 



TABLE III: Stability of the IS fit result, cr^^eory i^^fers to whether or not the theory error is included 
in the fit. 
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FIG. 3: Fit results for and Ai on the left, and and \Vcb\ on the right. The fit to the 
B XcCy data only {B — > Xciv and B Xg^ data combined) is shown by a dashed blue line 
(solid red line). The regions correspond to Ax^ = 1. 

IV. KINETIC SCHEME ANALYSIS 
A. Theoretical input 

Spectral moments oi B ^ Xciiy decays have been derived up to (9(l/m3) in the kinetic 
scheme [8] . Compared to the original paper, the theoretical expressions used in the fit contain 
an improved calculation of the perturbative corrections to the lepton energy moments [16] 
and account for the -Emm dependence of the perturbative corrections to the hadronic mass 
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moments [17]. For the B — > Xs^ moments, the (biased) OPE prediction and the bias 
correction have been calculated [9] . 

All these expressions depend on the b- and c-quark masses mb{n) and mc{fj,), the non- 
perturbative parameters and fJ^cin) ((9(l/m^)), p%{fi) and p\si.lA (^(V^^b))) ^.nd 

as [18]. The theoretical uncertainties can be separated into two categories: non-perturbative 
(related to the expansion in l/rrib) and perturbative (related to the expansion in ag). 

Following the recipe in Ref . [8] , the non-perturbative uncertainties are evaluated by vary- 
ing ji^ and //^ (p|) and p\g) by ±20% (±30%) around their "nominal" values = 4.6 GeV, 
mc = 1.18 GeV, //^ = o.4 GeV^, p% = 0.1 GeV^ i^l = 0.35 GeV^ and pis = -0.15 GeV^ 
corresponding to the uncertainty of the respective Wilson coefficient. All these variations 
are considered uncorrclatcd for a given moment. The theoretical covariance matrix is then 
constructed by treating these errors as fully correlated for a given moment with different 
-E'min while they are treated as uncorrelated between moments of different order. 

For the moments of the photon energy spectrum, we take 30% of the absolute value of the 
bias correction as its uncertainty. This additional theoretical error is considered uncorrelated 
for moments with different -Emin and different order. 

The theoretical uncertainties mentioned so far (non-perturbative, bias correction) arc used 
to construct a theoretical covariance matrix for the measurements and are thus included in 
the fit. The perturbative uncertainties are estimated by repeating the fit, setting ag to 
a different value. For lepton and photon energy (hadronic mass) moments, we vary ag 
within ±0.04 (±0.1) around the central value 0.22 (0.3). These ranges of variation follow 
the recommendations in Ref. [8]. The different treatment of the hadron mass moments is 
due to the fact that the calculation of the perturbative corrections to these moments is less 
complete. 

B. The function 

We use a function with seven free parameters: the semileptonic b ^ c branching 
fraction B{B X^lv), m^, mc, pi, p^, Pa and pig, 

x' = - (^>f")cov,^'((^)r"- - (^)f ) ■ (9) 
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Here, (X)^^- are the measured moments. {X)^"^ are the corresponding kinetic scheme 
predictions that depend on these free parameters. The covariance matrix is the sum of the 
experimental and theoretical error matrices. 

We determine the CKM element \Vcb\ by treating it as an eighth free parameter in the 
fit. iV^cbl is related to the semileptonic width r(i? — >• Xcii') by [5] 

T/ I / 1 P^Pi \ "'"/^ 

I r{B ^ XJu)-^-^] X (1 - 0.0018) X (1 + 0.30(q;,- 0.22)) (10) 



0.0417 V 0-105 , 

x(l - 0.66(m6 - 4.6 GeV) + 0.39(mc - 1.15 GeV) 

+0M3{i4 - 0.4 GeV^) + Om(pl - 0.1 GeV^) 

+0.05(//^ - 0.35 GeV^) - 0.01 (pi^ + 0.15 GeV=^) . 



Using this expression, we calculate r(S — > Xciu) from \Vcb\ and add the following term to 
the function, 

X -X +( riB^xju) -^'^ 

As fiQ and p^g are determined from B* — B mass splitting and heavy quark sum rules 
and because the expressions depend only weakly on these parameters, we fix Pq and p\g to 
0.35 ± 0.07 GeV^ and -0.15 ± 0.1 GeV^, respectively, by adding the following terms to the 

function, 

^12 ^ ^2 ^ ^^2^ _ Q_35 GeV2)V(0.07 GeV')^ + {pl^ + 0.15 GeV3)V(0.1 GeV^)^ . (12) 
The minimization of the X"^ is performed using MINUIT [15]. 



C. Fit results and discussion 



The result of the kinetic scheme analysis is shown in Table IV and in Figs. 4 and 5. The 
value of the x^ function at the minimum is 17.76, compared to (31 — 7) degrees of freedom. 
All results are preliminary. 

To assess the stabihty of the fit, we have repeated the analysis using lepton energy 
moments only, hadron mass moments only and photon energy moments only. The result 
is shown in Fig. 6. In general, changes are well covered by the fit uncertainty though the 
B — > data seems to prefer lower values of m^. 

Finally, the result for \Vcb\ reads 

iKbl = (41.93 ± 0.65(fit) ± 0.07(a,) ± 0.63(th)) x 10"^ . 
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Bxj. (%) 


rrib (GeV) m 


•c (GeV) 


pI (GeV2) 


pI (GeV3) pI (GeV2) 


Pis (GeV3) 


value 


10.590 


4.564 


1.105 


0.557 


0.162 


0.358 


-0.174 


(7(fit) 


0.164 


0.076 


0.116 


0.091 


0.053 


0.060 


0.098 


a{as) 


0.006 


0.003 


0.005 


0.013 


0.008 


0.003 


0.003 




1.000 


-0.023 


0.003 


0.229 


0.192 


-0.147 


-0.024 


mb 




1.000 


0.983 


-0.729 


-0.623 


-0.024 


-0.111 


rric 






1.000 


-0.716 


-0.633 


-0.124 


-0.033 


pI 








1.000 


0.851 


0.005 


0.052 












1.000 


-0.046 


-0.156 


P-G 












1.000 


-0.071 


Pls 














1.000 



TABLE IV: Results of the kinetic scheme fit. The error from the fit contains the uncertainties 
related to the experiment, the non-perturbative corrections and the bias correction. a{as) is the 
uncertainty related to the perturbative corrections. In the lower part of the table, the correlation 
matrix of the parameters is shown. 

The first error is due to all uncertainties taken into account in the fit (experimental error in 
the moment measurements, non-perturbative corrections and bias correction to the moments, 
uncertainty in tb). The second error is obtained by varying as in the expressions for the 
moments and in Eq. 10. In Eq. 10, we vary by ±0.008 around the central value of 
0.22 [19]. The last error is a 1.5% uncertainty due to the hmited accuracy of the theoretical 
expression for the semileptonic width, assessed in Ref. [5]. 

V. SUMMARY 

We have performed a fit to the Belle measured spectral moments of the lepton energy 
and hadronic mass spectrum in charmed semileptonic B decays, and the photon energy 
spectrum of inclusive radiative B decays using expressions for the moments in terms of 
HQE parameters in the 15" mass and kinetic mass schemes. The fits produce values of \Vcb\ 
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Belle preliminary 




FIG. 4: Partial branching fractions and lepton energy moments, compared to the kinetic scheme 
fit result. The yellow bands show the theoretical uncertainty included in the fit (non-perturbative 
corrections, bias correction). The open symbols correspond to measurements not used in the fit. 

that are consistent between the two schemes. In the IS scheme analysis we find \Vcb\ = 
(41.49 ± 0.52(fit) ± 0.20(rB)) x 10-^ and in the kinetic scheme we obtain = (41.93 ± 
0.65(fit) ± 0.07(a;s) ± 0.63(th)) x 10^'^. The heavy quark parameters, m^™'^'^ and Ai, /i^, 
have been extracted with values that are consistent with previous determinations [2, 3]. 
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Belle preliminary 




FIG. 5: Same as Fig. 4 for the hadron mass and photon energy moments. 

Constant feedback between theory and experiment should further confirm the understanding 
of the OPE in all measurable regions of phase space. The accuracy achieved by the Belle 
measurements is unprecedented by any other experiment and the uncertainty on the heavy 
quark parameters and \Vcb\ reflect these improvements. 
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FIG. 6: Stability of the kinetic scheme fit. The fit is repeated using lepton energy moments only, 
hadron mass moments only and photon energy moments only. The ellipses are Ax^ = 1. 
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